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Incremental  Diffraction  Coefficients  for  the  Truncated  Half-Plane 
and  the  Calculation  of  the  Bistatic  Radar  Cross  Section  of  the  Disk 

1.  INTRODUCTION 

Exact  expressions  for  incremental  diffraction  coefficients  (IDC's),  also  called  equivalent  edge 
currents,  have  been  derived  for  the  perfectly  conducting  wedge  by  integrating  the  wedge  currents  over 
an  Increment  of  the  wedge1-2  and  more  recently,  for  any  planar  current  surface  or  combination  of 
planar  current  surfaces  by  direct  substitution  of  the  two-dimensional  (2-D)  far  fields  of  each  planar 
surface  into  convenient  general  expressions.3-4,5  Incremental  diffraction  coefficients  obtained  for  the 
total  current,  the  physical  optics  current,  or  the  nonuniform  current,  that  is,  the  difference  between 
the  total  and  physical  optics  currents,  are  referred  to  as  the  GTD  (geometrical  theory  of  diffraction),  PO 


(Received  for  Publication  14  November  1989) 
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(physical  optics)  and  PTD  (physical  theory  of  diffraction)  incremental  diffraction  coefficients, 
respectively.3  For  the  perfectly  conducting  wedge,  Mitzner1  concentrated  on  the  PTD  incremental 
diffraction  coefficients,  Michael!  derived  the  GTD  coefficients  initially2  and  later  the  PTD 
coefficients6,  and  Knott7  determined  the  PO  incremental  diffraction  coefficients  from  the  difference 
between  the  GTD  and  PTD  coefficients. 

Assume  incremental  diffraction  coefficients  have  been  determined  for  a  canonical  2-D  scatterer 
of  finite  cross  section.  The  integration  of  ihest.  finite  inciemental  diffraction  coefficients  (multiplied 
by  the  Incident  fle’d)  along  the  bounding  curves  of  three-dimensional  (3-D)  scatterers  produces  a 
uniform  high-frequency  solution  even  at .  ,i  near  caustics  where  conventional  stationary-phase 
high-frequency  techniques  fall.  As  an  example.  Shore  and  Yaghjian  obtained  uniform  IDC's  for  the 
narrow  (electrically  small)  strip  and  slit  and  used  them  to  calculate  the  effect  of  cracks  between  the 
panels  of  reflector  antennas  on  tin  antennas'  co-pularized  and  cross-polarized  far  fields  3,6 

For  canonical  problems  of  Infinite  cross  section,  like  the  wedge,  the  incremental  diffracMcn 
coefficients  contain  the  same  discontinuities  and  singularities  as  the  far-field  of  the  canonical  2-D 
problems  Specifically,  the  GTD  incremental  diffraction  coefficients  of  the  wedge  contain  the 
discontinuities  in  the  diffracted  H-wave  fields  at  the  face  angles  of  the  wedge,  the  singularities  in  the 
E-  and  H-wavc  far  fields  at  a  face  angle  of  the  wedge  when  an  incident  plane  wave  giazes  that  face  from 
the  outside,  and  the  singularities  in  the  total  diffracted  far  fields  at  the  shadow  and  reflection 
boundaries  Integrating  the  PO  current  separately  and  using  PTD  incremental  diffraction  coefficients, 
instead  of  GTD  incremental  coefficients,  eliminates  the  singularities  at  the  shadow  and  reflection 
boundaries,  but  retains  the  face-angle  discontinuities  and  singularities.6 

Figure  1  shows  the  specularly  scattered  fields  of  a  perfectly  conducting  infinitely  thin  disk 
computed  by  numerically  Integra* ing  the  PTD  incremental  diffraction  coefficients  for  the  infinite 
half-plane  arcund  the  rim  of  the  disk,  and  adding  the  result  to  the  PO  fields.  The  comparison  in 
Figure  1  of  the  PO  plus  ID*"1  fields  with  the  PO  fields  alone,  and  with  the  exact  scattered  fields  obtained 
from  the  eigenfunction  solution  to  the  disk,  shows  that  the  in'egrated  FID  incremental  coefficients  of 
the  hall- plane  become  divergent  at  grazing.  (Interestingly,  for  specular  scattering  from  the  disk,  the 
PO  and  integrated  PTD-IDC  far  fields  are  identical  for  the  two  incidental  plane-wave  polarizations.) 
Thus,  Figure  1  illustrates  the  limitations  of  infinite  half-plane  IDC’s  and  motivates  tills  report’s 
primary  objective:  to  determine  accurate  high-frequency  incremental  diffraction  coefficients  for  the 
leading  and  trailing  edges  of  the  truncated  half-plane  in  order  to  eliminate  the  singularities  and 
discontinuities  in  the  diffraction  coefficients  of  the  Infinite  half-plane.  Of  course,  the  2  D  truncated 
half- plane,  that  is  the  strip,  has  an  exact  eigenfunedon  solution  that  could  be  used  to  find  < he  3-D 
incremental  coefficients  by  substituting  the  exact  solution  into  the  general  expressions  of  Reference  3. 
However,  the  summation  of  eigenfunctions  takes  a  considerable  amount  of  computer  time  for 


6.  Michaeli,  A.  (1386)  Elimination  of  infinities  in  equivalent  edge  currents,  Part  I:  fringe  current 

components,  IEEE  Tram,.  Antennas  Fropagai.  AP-34:yt2-yi8. 

7.  Knott,  E.F.  (1985)  The  relationship  between  Mitzner’s  ILDC  and  Mlchaeli’s  equivalent  currents, 

IEEE  Trans.  Antennas  Propagat.  AP-33: 112-114. 
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Figure  lb.  Specular  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  for  Plane- Wave  Illumination  with  the  Electric  Field 
Polarized  Parallel  to  the  Plane  of  Incidence;  - :  Exact, - :  PO, - —  —  •-  :  PO  +  Half-Plane  IDC 


electrically  large  strips  and  the  available  computer  codes  encountered  numerical  difficulties 
evaluating  the  required  Mathieu  functions  when  the  strips  became  wider  than  about  5  wavelengths. 
Moreover,  the  exact  strip  solution  yields  the  incremental  diffraction  coefficients  for  rectangular 
incremental  strips  only,  rather  than  for  general  trapezoidal  strips  that  are  able  to  conform  to  the  rim 
of  an  arbitrarily  shaped  polygonal  plate  or  disk.  Thus,  we  found  it  desirable  to  determine  an  accurate, 
uniform  high-frequency  solution  to  the  perfectly  conducting  truncated  half-plane  with  the  leading 
and  trailing  edges  treated  separately  (though  not  necessarily  independently)  so  that  IDC's  for 
trapezoidal  strips  could  be  found. 

Previous  IDC  work  has  also  involved  high-frequency  solutions  to  the  strip.8*9-10  Coleman  et  al8 
combined  the  use  of  rectangular  incremental  strips  that  do  not  conform  to  the  rim  of  the  general  3-D 
plates  to  which  they  are  applied,  with  an  early,  relatively  inaccurate,  high-frequency  2-D  strip 
solution  of  Ufimtsev. 1 1  Sikta  et  al9  also  used  nonconformal  rectangular  incremental  strips  and 
obtained  scattering  in  the  plane  of  incidence  only.  Michacli10  used  conformal  incremental  strips,  but 
evaluated  the  current  integrals  of  the  incremental  strips  asymptotically  so  that  singularities  at 
grazing  incidence  arise.  None  of  these  previous  strip  solutions  included  diffraction  of  the  nonuniform 
leading  edge  currents  at  the  trailing  edge  of  the  strip.  An  essential  feature  of  the  high-frequency 
solution  derived  in  this  report  is  that  the  diffraction  at  the  trailing  edge  of  the  H-wave  nonuniform 
leading  edge  current,  which  becomes  appreciable  near  grazing  incidence,  is  taken  into  account  even 
when  the  leading  and  trailing  edges  of  the  incremental  strip  are  not  parallel. 

The  analysis  begins  by  finding  convenient  expressions  for  two-dimensional  fields  radiated  by 
the  truncated  leading  and  trailing  edge  currents  of  a  half-plane  illuminated  by  a  TM  or  TE  plane  wave. 
These  2-D  leading  and  trailing  edge  diffracted  fields  are  then  used  to  compute  the  high-frequency 
bistatic  radar  cross  section  of  an  infinite  strip.  This  high-frequency  solution  is  then  compared  with 
the  exact  eigenfunction  solution  for  the  strip. 

Next,  the  2-D  diffracted  fields  of  the  truncated  half-plane  are  inserted  into  the  general 
expressions  of  Reference  3  to  obtain  the  PTD  incremental  diffraction  coefficients  for  the  leading  and 
trailing  edges  of  the  truncated  half-plane.  Finally,  these  PTD  incremental  diffraction  coefficients  are 
applied  to  finding  the  approximate  far  fields  scattered  by  a  circular  disk  illuminated  with  a  TM  or  TE 
plane  wave.  These  approximate  high-frequency  scattered  fields  are  compared  with  the  exact  far  fields 
computed  from  the  eigenfunction  solution  to  the  disk. 


8.  Coleman,  J.R.  (1973)  Investigation  of  Radar  Scattering  by  Simulated  Aircraft  Duct/ Engine 
Combinations,  Air  Force  Avionics  Laboratory,  AFAL-TR-73-361. 

9  Sikta,  F.A.,  Burnside,  W.D.,  Chu,  T.T.,  and  Peters,  L.,  Jr.  (1983)  First-order  equivalent  current 
and  comer  diffraction  scattering  from  flat  plate  stuctures,  IEEE  Trans.  Antennas  Propagat. 
AP-31:584-589. 

10.  Mlehaeli,  A.  (1987)  Equivalent  currents  for  second-order  diffraction  by  the  edges  of  perfectly 

conducting  polygonal  surfaces.  IEEE  'Dans.  Antennas  Propagat.  AF-35: 183-190. 

1 1.  Ufimtsev,  P.  Ya.  (1958)  Secondary  diffraction  of  electromagnetic  waves  at  a  strip.  Journal  of 

Technical  Physics  28  (No.  3). 
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2.  THE  FAR  FIELDS  RADIATED  BY  THE  TRUNCATED  CURRENTS  OF  A 
HALF-PLANE  ILLUMINATED  BY  A  PLANE  WAVE 

The  general  expressions  in  Reference  3  yield  the  3-D  incremental  diffraction  coefficients  of  a 
planar,  perfectly  conducting  2-D  canonical  scatterer  through  direct  substitution  of  the  2-D  far  fields  of 
this  canonical  scatterer.  Thus  we  begin  the  derivation  of  IDC's  for  the  truncated  half- plane  by 
deriving  the  far  fields  radiated  by  the  currents  on  a  truncated  half-plane,  or  more  precisely  the  far 
fields  of  the  truncated  currents  of  an  infinite  half-plane  under  plane-wave  illumination.  Figure  2 
shows  the  truncated  half-plane  lying  between  x  =  0  and  x  =  s  and  stretching  from  to  in  the 
/.-direction.  The  spherical  angles  (<J>,0)  of  the  observation  direction  r,  and  the  spherical  angles  (^.OJ  of 
the  propagation  direction  kof  the  incident  plane  wave  are  also  shown  in  Figure  2. 

OBSERVATION  OBSERVATION 


Figure  2.  Geometry  of  the  Truncated  Half-Plane  Illumiiicu.- A.  by  a  Plane  Wave 


The  far  fields  radiated  by  a  current  K  on  the  truncated  half-plane  are  given  by  the  integral^ 

s 

■772  ikr0  x  J  K(x')e-Ikx' sln  °°  cos<’>  dx'  (la) 


H(r ) 


r  -><*> 


cikr0  •  fglTt/4 
(8xkp  sin  0O) ! 


E(r )  — Z0r0  x  I-I(r ),  r ->oo  ,  (lb) 

where  the  unit  vector  r0  in  Eq.  (1)  is  defined  as  the  unit  vector  r  evaluated  at  <{)  and  0  =  n  -  0O,  that  is, 
r0  =  p  sin  0O  -  /.  cos  0O.  The  exp(-icot)  time  dependence  has  been  suppressed  and  Z0  is  the  plane-wave 
impedance  of  free  space.  For  a  normally  incident  plane  wave  (0O  =  90°),  Eq.  (1)  reduces  to 
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E(p)  — Z0px  H(p) ,  p  ->oo  .  (2b) 

The  far-field  expressions  In  Eq.  (2)  for  normal  incidence  are  especially  important  because  once  they 
are  found,  the  far  fields  for  an  obliquely  incident  TM  plane  wave  are  determined  immediately  by 
substituting  k  sin  0O  for  k  and  multiplying  by  exp(-lkz  cos  80).  The  H0  (or  E^)  far-field  for  an  obliquely 
incident  TE  plane  wave  is  also  determined  by  this  simple  substitution  and  multiplication.  However, 
the  (or  E0)  far  field,  which  occurs  only  at  oblique  incidence  in  the  TE  case,  must  be  obtained  by  an 
alternative  procedure;  for  example,  by  integrating  the  current  in  Eq.  (I).3 


2.1  The  TM  (Parallel  Polarized)  Fields:  Normal  Incidence 

The  total  current  on  a  perfectly  electrically  conducting  infinite  half-plane  illuminated  by  a 
normally  incident  (0O  =  90°)  TM  plane  wave,  given  by 

E,  =  zE|  exp  [ikp  cos  (<J>  -  4>0) ]  .  (3) 

may  be  found  from  Eq.  (8. 14)  of  Bowman,  Senior  and  Uslenghi12: 

K  =  zKz  =  z^E t  sin  tjp  ei7t/4 

.  (^J/2e‘kx  -2iV2~  cos  ^  e-lkxcos^F(^2kcos2^xj  (4) 

where  F  is  the  Fresnel  integral  defined  by 
x 

1  re11 

F(x)  =  -=  -=  dt  .  (5) 

V2x  J  Vt 

0 


Substitute  the  to<al  current  K  from  Eq.  (4)  into  the  integral  of  Eq.  (2),  and  make  use  of  the  following 
integral  of  the  Fresnei  integral  derived  in  Appendix  A, 


x 


(F°| a  +  1  lx)  ' 
|  a  +  1 1 1/2 


(6) 


12.  Bowman,  J.J.,  Senior,  T.B.A.,  and  Uslenghi,  P.L.E.  (1969)  Electromagnetic  and  Acoustic 
Scattering  by  Simple  Shapes,  Amsterdam:  North-Holland. 
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where  F°(U  denotes  F(t)  when  t  >  0  and  the  complex  conjugate  of  F(t)  when  t  <  0,  to  obtain  the  TM  far 
fields  of  the  total  current  for  normal  incidence 


E™(p)P^T  zE, 


A  elkP  sm  2 


(xkp) l/2  cos  <j>  0  +  cos  4> 


V2(l  -  cos  ) i .  . 

•  ...  ==r  F°(ksl  1  -  cos  <{>  I ) 

1 1  -  cos  4>  | 

-2  cos^  e-,ks(cos<>o  +  c03^F^2kscos2^  , 


0  <  4>.  4>0  <  2jc  , 

H™(p)  ~  p  x  E™(p)  .  p->°°  .  (7b: 

(We  retain  the  absolute  value  signs  in  1 1  -  cos  4>|  because  when  the  IDC’s  are  found  by  substitution  of 
Eq.  (7)  into  the  general  expression  of  Reference  3,  cos  <J>  is  replaced  by  a  variable  that  may  become 
greater  than  one;  see  Section  3.) 

Eqs.  (7)  give  the  far  fields  of  the  total  current  emanating  from  either  a  leading  or  trailing  edge  of 
the  truncated  half-plane.  To  obtain  the  TM  far  fields  of  the  truncated  nonuniform  current,  simply 
subtract  the  TM,  PO  far  fields, 


ErofrM»(p)P^aEl  ^2 


.  ‘j’o  ‘t’o 

;i(kp  +  jr/4)  s*n  2  C0S  2 
(jikp)  cos  <(>0  +  cos  4 


_ [  T  ,  iks(cos  4> o  +  cos  <J>)1 

0  +  COS  4)  L  J 


H^'tpl^pxE^/Zo. 


from  Eq.  (7)  to  get 


T  0 

En"m,,(5,p~  ?-E-  cosvT-^oVc  F0(ksl1  - C0SI>|1 

-  I  cos  ^  |  [2e-‘ ks(cos  +  cos  w  F^2ks  coS2  ^ 

+  V2"e1't/4  ( 1  -  e■iks(cos,:,0  +  cos<f>))  j, 


0<4>.  4>o < 


H  nufrM)(p)  ~  p  x  E  nu,™’(p)/Z0  ,  p->oo  . 
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Depending  on  whether  the  incident  plane  wave  makes  an  angle  4>0  greater  than  or  less  than 
90"  (or  less  than  or  greater  than  270°),  Eq.  (9)  gives  the  far  fields  of  the  nonuniform  half-plane  current 
of  the  leading  or  trailing  edge,  respectively,  truncated  at  the  trailing  or  leading  edge,  respectively.  Of 
course,  In  reality,  there  will  exist  multiple  interactions  between  the  leading  and  trailing  edges 
dominated  by  the  nonuniform  current  from  the  leading  edge  diffracting  at  the  trailing  edge.  However, 
for  the  TM  Incident  fields  the  nonuniform  current  decays  as  the  inverse  of  the  square  root  of  the 
distance  from  the  edge,  and  thus  Eq.  (9)  is  a  good  approximation  for  the  nonuniform  far  fields  of  the 
trailing  as  well  as  the  leading  edges  of  electrically  large  strips  (s  >  X). 


2.2  The  TM  Fields:  Oblique  Incidence 


As  mentioned  in  Section  2,  the  far  fields  radiated  by  the  total,  PO,  and  nonuniform  currents  can 
be  determined  for  an  obliquely  incident  TM  plane  wave  essentially  by  replacing  k  with  k  sin  0o  in  the 
far  fields  for  normal  incidence  and  multiplying  by  exp(-ikz  cos  Go).  In  particular.  Eqs.  (7)  generalize 
for  oblique  incidence  to 


E  nufrM)(r ) 


-SgEi 


C0  2^2 


-in/4 


sin 


$0 


ik 


COS  ({>0  +  cos  <j> 


V2~(l  -  cos  <t>) 
V 1 1  -  cos  4>  | 


F°(ks  sin  00 1 1  —  cos  4>  I ) 


2e~i  ks  sin  90  (cos  +  cos  $) 


F^2ks  sin  0O  cos2 


+  V2’ei,t^4 


-1  ks  sin  G0  (cos  +  cos 


(10a) 


0  <  <  2n, 


0  £  0O  <  k  , 

H  nu(TM,(r )  ~  r  0  x  E  nufIM)(r)/Z0,  r->~  (10b) 

where  0§  is  the  unit  vector  0  evaluated  at  0  =  n  -  0O,  and 

eik(p  sin  O0  -  z  CCS  e0)  gin/4  ^ 

^°-  (8xkp  sin  0q) J/2  ' 

In  Section  3,  the  2-D  TM  fields  [Eqs.  (10)1  are  substituted  into  Eqs.  (23)  and  (24)  of  Reference  3  to  obtain 
the  3  D  incremental  uiffi action  coefficients  for  TM  plane  wave  incidence  on  the  leading  and  trailing 
edges  of  the  truncated  half-plane. 
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2.3  The  TE  (Perpendicular  Polarized)  Fields:  Normal  Incidence 

The  2-D  far  fields  of  a  perfectly  conducting  truncated  half-plane  for  a  normally  incident  (60  =  90°) 
TE  plane  wave,  given  by 

H,  =  zHj  exp[-ikp  cos(4>  -  <i>o51  U2) 

begins,  as  in  the  TM  case,  with  the  total  current  on  the  infinite  half-plane: 

K  -  Kxx  =  x  H,  sign  (re  -  4>0)2V2~e"i  3/4  e'Ik  cos  F|2k  cos2  xj  (13) 

which  can  be  found  from  Eq.  (8.31)  of  Bowman,  Senior  and  Uslenghi.12  The  sign  (re  -  <t>0)  function  in 
Eq.  (13)  is  +1  or  -1  for  4>0  less  than  or  greater  than  re,  respectively.  Substituting  K  from  Eq.  (13)  into 
Eq.  (2),  and  making  use  of  Eq.  (6)  to  perform  the  Integration  of  the  Fresnel  integral,  yields  the  TE  far 
fields  of  the  total  current  for  normal  incidence 


Hsra(p)  PZT 


sin  <b 


TT  I  X  \  V'  W 

/j  1  s  75  ~  <*>0  (rckp) cos  <J>o  +  cos  4* 


.  V2 

-Iks  (cos  $0  +  cos<>)  ( olM,  *-£.2  22 1  _ _ 

4>0 

COS-g- 

_  rnM.*  !  i  a  1  \ 

e  ,  to  j,  ^2kS  COS'4  -jjf j  - 

0  S  (j>,^0  <  2re, 


EsTC(p)  ~  -Z0  P  x  HsTC(p),  p->~  (14b) 

(Again,  we  retain  the  absolute  value  signs  in  1 1  -  cos  <J>!  because  when  the  IDC’s  are  found  by 
substitution  of  Eq.  (14)  into  the  general  expressions  of  Reference  3,  cos  <&  is  replaced  by  a  variable  that 
may  become  greater  than  one:  see  Section  3  below?.) 

The  TE  far  fields  of  the  truncated  nenuniform  current  are  determined  by  subtracting  the  TE,  PO 
far  fields. 


_  PO(iti)  _  a  .  .  e1(kp+ ^  sin  4. 

H  (P)  P2-  -  ZH,  Slgnu  -  ♦„)  l2lkp),/2  C03  4>0+  cos* 

f  j  _e~i  ksfeos  <p0  +cos  4>)  1 


(15a) 


_  TOflB)  _  TOfTE)  _ 

E  (p)  — Z0  pxH  (p)  ,  p->°° 


from  Eq.  (M)  to  get 


(15b) 


10 


_  inifTE)  a  sin  6 

H  (p)  zH,  sign(?t  -  4> 0)  7—5 — 7772  - 1 - T 

K  ~  1  b  0  (rtkp)1/2  cos  ij>0  +  cos  <J> 


V2  I  COS  7? 


e-i  ks  (cos  0O  +  cos  0)  p^2ks  cos2  &  j  _  _ L— -  -  jL  F°(ks  I  1  -  COS  4>  I )  (16a) 


elr‘/4 

44 


1  -  e'1 


-1  ks  (cos  <J>0  +  cos  <1 


’] 


f  .  0  2  <p.  4)0<  2tc, 


_  nu(TE) 


_  nu(TE) 

(p)  ~  -Z0  p  x  H  (p)  ,  p~>«  . 


(16b) 


Eqs  (16)  present  a  good  high  frequency  approximation  to  the  fields  radiated  by  truncated  TE 
nonuniform  current  emanating  from  the  leading  edge  for  all  angles  of  incidence.  However,  unlike  the 
analogous  TM  Eqs.  (9),  Eqs.  (16)  need  t )  be  modified  for  the  trailing  edge,  when  the  incident  angle  is 
near  grazing  (<+>0  =  0°),  to  take  into  account  the  nonuniform  current  that  emanates  from  the  leading  edge 
and  diffracts  at  the  trailing  edge.  Fortunately,  this  modification  required  for  the  trailing  edge  under 
TE  illumination  is  accomplished  through  a  simple  multiplicative  factor  that  will  be  derived  in 
Section  2.5. 

2.4  The  TE  Fields:  Oblique  Incidence 

As  mentioned  in  Section  2,  the  H0  or  far  fields  radiated  by  the  nonuniform  current  can  be 
determined  for  an  obliquely  incident  TE  plane  wave  essentially  by  replacing  k  with  k  sin  0O  in  the  far 
fields  for  normal  incidence  and  multiplying  by  exp(-ikz  cos  0O).  In  particular,  Eqs.  (14)  and  (16) 
generalize  respectively  for  oblique  incidence  to 


Hfso’(r)  0gr2“  -  OqH,  sign(7t  -  4>0) 


C0242g 


-t  Tt/4 


sin  1)) 


Ik 


cos  <t>0  +  cos  ‘t1 


a-l ks  sin  90(cos 


0+  cos  <»  F^2ks  sin  0O  coS2 


-42 


$0 


cos  ~2 


jl“  cosTp2'  F°(kSSln  0o 1 1  -  cos  4»  I ) 


(17a) 


0  <  ,4>0  <  27C , 

0  <  0O<  7t, 

El(r)$  r72”  -Z0  r  0  x  H~q  (r)^g,r->«. 


(17b) 


and 
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■Vw.  Mr -if  \f .W  HA  iisVjes* 


HonUnE)(r)  og  -  og  H,  signtTC  -  <j»0)  ^  co“®‘"+^ 


(|)q  +  COS  ({> 


i  V2"e-,II/4 


e-iks  sin  O0(cos«0  +cos  4)  p^ks  sin  Qq  cos2 


V2 


COS 


4*0 


1  -  COS  <J)  I 


Yj2  F°(ks  sin  00 1 1  -  cos  <j>  I ) 


j  -1  ks  sin  0o(cos  t/o  +  cos  $) 


(18a) 


0  <  4> .4>o  <  2k, 


_  mini) 
E « 


(r) 


$ 


A  _  nj(TE) 

i  rox  H  o  (r) 


A 


0 


jt 

0  • 


r->o°  . 


0  <  0O  <  71, 


(18b) 


For  the  TE  case,  there  may  also  exist  H0  and  E0  far  fields  when  the  plane  wave  is  obliquely 
incident.  These  H^,  and  E0  far  fields  can  be  found  by  integrating  the  TE  current  that  exists  on  the  half¬ 
plane  for  oblique  incidence.  In  particular,  the  nonuniform  far  fields,  H0nuri  E)  or  E0nu(TE*,  can  be  found 
for  the  truncated  half-plane  by  subtracting  the  PO  far  fields,  H0POtTE)  or  E0POflE),  of  the  truncated  half- 
plane  from  the  TE  far  fields,  Hs<;i  or  Es0  ,  obtained  by  integrating  the  total  current  over  the 
truncated  half-plane.  Before  we  do  these  integrations  to  determine  the  or  Ee  fields,  however,  let  us 
address  an  objection  that  may  arise  at  this  point  in  the  analysis. 

If  the  total  TE  current  that  we  are  integrating  were  the  exact  total  current  of  the  2-D  perfectly 
conducting  strip,  then  the  basic  theory  of  2-D  scattering  from  perfect  conductors  (see  Chapter  1  of 
Reference  12),  tells  us  that  these  total  currents,  when  integrated  will  not  produce  an  H0  and  E0 
component  in  the  far  fields.  Of  course,  we  are  not  integrating  the  total  currents  of  the  strip,  but  the 
truncated  currents  of  the  Infinite  half-plane.  Thus,  we  cannot  expect  the  H6  or  E0  far  fields  to  be  zero. 
Still,  one  could  argue  that  when  we  eventually  apply  these  far  fields  to  both  the  leading  and  trailing 
edges  of  Incremental  strips,  the  combined  truncated  currents  from  both  the  leading  and  trailing  edges 
should  closely  approximate  the  exact  total  currents  of  the  strip,  and  thus,  we  should  ignore  the  and 
E„  components  ofTF  far  fields  at  all  points  in  the  analysis.  This  argument  is  valid  provided  the 
leading  and  trailing  edges  of  the  lncr  mental  strip  arc  parallel.  Only  for  these  icclangular  stilps  is  the 
incremental  strip  part  of  an  infinite  2-D  strip  for  which  the  above-mentioned  theory  in  Chapter  1  of 
Reference  12  holds.  For  trapezoidal  incremental  strips,  that  is,  strips  with  nonparallel  leading  and 
trailing  edges,  we  must  include,  in  general,  the  TE  and  E0  far  fields  radiated  by  the  truncated  half- 
plane  currents  for  oblique  Incidence. 

To  obtain  the  total  scattered  TE  or  E0  far  fields  it  suffices  to  obtain  the  z-component  of  the 
total  current  on  the  half- plane  since  then  Eqs.  (9b)  and  (14b)  of  Reference  3  combined  with  Eq.  (17a) 
already  derived  for  IIs9  yield  HS(J).  Since  Kz  on  the  surfaces  of  the  half-plane  is  equal  to  Hx  within  a 
sign,  we  seek  an  expression  for  Hx  on  the  half-plane  surfaces. 
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For  the  normally  incident  TE  plane  wave  given  bv  Eq.  (1 1).  Eq.  (8.28)  of  Reference  12  gives  the 
total  rl7  everywhere  in  space.  A  general  procedure  for  generating  the  solution  for  oblique  incidence  is 
given  in  Reference  12.  For  H-polarization  lot  V  be  obtained  from  the  normal  incidence  (two- 
dimenslona1;  solution  by  substituting  h  ->  k  sin  80  and  multiplying  by  exp(-ikz  cos  0O).  Then  Eq.  (1.6) 
of  Reference  12  adapted  to  the  form  of  the  obliquely  incident  plane  wave 

Hinc  =  -H,  (cos  <i>0  cos  G0  x  +  sin  4>0  cos  60  y  -  sin  0O  z) 


•  exp  -ik  (x  cos  <{>0  sin  0O  +  y  sin  <j>0  sin  0O  +  z  cos  0o)  j 


gives  the  total  magnetic  field 

_  i  cos  0n  /3V  a  9V  an  a 

H=-k^fcx+3Fyysinoovz 


with,  from  Eq.  (8.28)  of  Reference  12, 


-1  Jt/4 

xr  -  u  - _  -i  kz  cos  O 

V-H'  ^  e 


°  jexp 


ik  sin  0op  cos  (<j> 


'2k  sin  0O  p  cos 


^  (4>  - 


exp  j^-ik  sin  0O  p  cos  (<J>  +  $0)  j  ^2k  sin  0op  cos  ^-(<J)  +  (J>0)  j 


y  ((o)  =  J  elt2  dt  . 


Then  Hx  on  the  surfaces  of  the  half-plane  is  given  by 

JL  cos  e0  3V 
x  "  k  sin  0O  dx 


V  =  2H,  e-1  kz  cos  00  exp(-ik  sin  0O  cos  (|)0x)  N/ik'shT^lTcos  (2 

where  the  -  or  +  sign  applies  on  the  top  and  bottom  half-plane  faces  respectively.  The  current  on  the 
top  and  bottom  faces  of  the  half-plane  is  then  given  by 

K(b)  =  +  Hv 


=  ± cos  ^ ccs  <p0  V  +  H,e 1,1 74 e'lkz 005  0  0  cos  0ocoslf- 


e'hsln0°x 
2  x'/~  • 


The  total  current  Is  obtained  by  adding  the  currents  on  the  two  faces: 
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Ki{  =  sign ( 7i  -  <j>0)H,  2V2"cos  O0e",kzoos  °° 

•  [e-u/4  e~lk  sln  °o  003 <>ox  cos  <J)0  F  ^2k  sin  0ocos2  ^x) 

] 


,ijr/4 


r - - — - -  —  V.UO  cy 

ynk  sin  0O  I 


,1k  sin  OqX 


(26) 


0  cos  <}>0  +  cos  6 


with  F(x)  given  by  Eq.  (5).  Integrating  K,101  over  the  half-plane  from  0  to  s  and  making  use  of  Eq.  (6) 
yields  [see  Eq.  (9b),  Reference  3). 

s 

Aqj,  =  Jiyot  (x)  e-lkx  sln0°oos*  dx 
0 

•  [cos  *0  e"lks  sln  “o'""*®*  “*>  f[21cs  sin  O0cos2  ^ 

I  I 

+  ^  I  1  -  COS  <)>  j  F°  (kS  Sin  °0  I  1  -  COS  $  I J  . 

From  Eq.  (14b)  of  Reference  3  we  then  obtain 

=  -  sin  O0  (C0A0k  +  cot  0O  cot  4>  Hs0) 
from  which,  with  Eq.  (17a), 

H™  (?)  $  r£“  sign { 7i  -  4»0)Ht  2%r2e-i*l*  yg  cos  O0 

.  e-l  ks  sin  O0(cos  ^0+  a»  0)  p  Sin  0O  C0S2  ^ . 

OS({i,<i)O<27r,OS0o<7i:, 

ES  (?)  V  -r0x  II™  (?)  $  •  r-><»  .  (29b) 

The  corresponding  PO  fields  are  readily  obtained  from  the  expressions  for  H^’0  for  the  faces  of 
the  infinite  wedge  given  on  pp.  32  33  of  the  technical  report  version  of  Reference  3  by  multiplying  by 
the  factor  j  1  -  expj^-iks  sin  0O  (cos  <)>0  +  cos  <)>)jj: 


(27) 


(28) 


(29a) 


PO(TE)  /  op  r  -i 

II  6  (?)  $  r^.“  -  sign(7i  -  <S>0)H,  -n~  cos  0OI  1  -  e_,ks  sin  °o(°os<!,c+  006  <•  |  . 


2C 

t  U 

l  XV 


0  <  ,(J»0  c  2ti,  0  <  0o<tc, 


(30a) 
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(30b) 


TOfTE)  a  a  POrTE)_  A 
Eo  (r)  0J~  -r0x  H?  (r)  $  .  r->°°  , 


and  hence,  subtracting  Eq.  (30)  from  Eq.  (29), 


H 


nu(TE) 


2Cr 


0  (r)  <J)  rzr  -  sign(x  -  <j> q) H (  cos  0o 


e~lks  sin  O0(cos  <>0+  cos  £) 


F  ^2ks  sin  0O  cos2 


♦[ 


1 


e-i  ks 


sin  O0(tos  <>0+  cos 


(31a) 


0  ^  4> ,  4>  c  ^  2rc ,  0  S  0O  <  x, 
nufrc)  A  A  nutTE)  A 

E0  (r)0g  —  r0x  Hd  (r)  <J»  .  r->«  .  (31b) 


2.5  Modification  of  the  TE  Fields  for  the  Trailing  Edge 

Consider  the  truncated  half- plane  in  Figure  2  with  a  TE  plane  wave  Incident  near  (j>0  =  180°.  The 
nonuniform  current,  defined  as  the  difference  between  the  total  and  PO  currents,  emanating  from  the 
leading  edge  decays  very  slowly  (unlike  In  the  TM  case)  toward  the  trailing  edge.  This  appreciable 
nonuniform  TE  current  impinging  upon  the  trailing  edge  produces  diffracted  fields  from  the  trailing 
edge  that  modify  the  diffracted  TE  fields,  Eq.  (18).  applied  to  the  trailing  edge.  We  can  determine  this 
modification  approximately  by  first  deriving  from  Eqs.  (20)  and  (21)  that  the  total  TE  current  for 
oblique  incidence  on  the  half-plane  is  given,  for  all  angles  of  incidence,  by 

K  =  -sign(7t  -  <}>0)2y  x  0oH,  e-'k (xsln  °o  003  «o+  z  cos  o0)  F^2kx  sin  0ocos2 

(j)  A 

A  ,  cfo”  COS  0O  TT 

+  2z  H,e1,t/4  —  2  elk(x  sln  0o  -  *  o0)  (32a) 

v  n  Vkx  sin  0O 

where  0o  is  the  unit  vector  0  at  0  =  0O.  For  2kx  sin  0o  »  1,  the  second  term  in  Eq.  (32a)  is  small 
compared  to  the  first  term,  and  Eq.  (32a)  can  be  written  approximately  as 

K  =  -sign(jt  -  4>0)2y  x  0oH,  e"lk(x  sln  °°  cos «0+  z  cos  e0) 

•  V2"e'1;t/4  F^2kx  sin  0ocos2  Tpl  .  2kx  sin  0O»  1  .  (32b) 
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The  TE  H-f)eld  corresponding  to  this  cuiTent,  Eq.  (32b),  impinging  on  the  trailing  edge  (x  =  s)  of  the 
truncated  half-plane  is 

HImp  =  -'Qq  h,  e‘ik(x  sln  °o  cos <)0+7.  n»o0)  .  2  \/2*e',n/4  F^2kx  sin  0ocos2  (33) 

2kx  sin  O0  >>  1  . 


Now,  the  obliquely  incident  TE  plane-wave  field  alone  along  the  half-plane  Is  given  by 

H,  =  -O0  H,  e"lk(x  sln  °°  cos<’o+z  cos0o)  (34) 

Comparing  Eqs.  (33)  and  (34),  we  see  that  the  Impinging  Held,  that  is  the  effective  incident  field  on  the 
trailing  edge  (x  =  s),  is  approximately  the  actual  incident  field  multiplied  by  the  factor 

2V2~e"llt/4  F^2ks  sin  0^  cos2  Ip  j,  (35) 

where  the  superscripts  L  on  <J>y  and  <J>0  in  (35)  emphasise  that  they  are  the  G0  and  <J>0  angles  defined  with 
respect  to  the  leading  edge  -  even  though  (35)  is  used  as  a  factor  multiplying  Eqs.  (18)  and  (31)  when  and 
only  when  Eqs.  (18)  and  (31)  are  applied  to  the  trailing  edge.  If  the  leading  edge  makes  an  angle  y  with 
the  incremental  strip,  so  that  s  is  given  in  terms  of  the  length  L  of  the  incremental  strip  by 

s  =  L  sin  y  ,  (36) 


(35)  becomes 


VsTe-11^4  Ff2kL  sin  y  sin  0q  cos2 


4*  o 


(37) 


If  in  addition  y  Is  taken  as  the  direction  of  propagation  in  the  half-plane  of  the  PO  current  or, 
equivalently,  of  the  total  and  nonuniform  current  near  grazing  (4>o  -  180°)  then  sln  y  equals 
-sin  Oq <J> o  yjl  ~  sin2  Oq  sin2  and  (37)  becomes 


V2Vl7t/4 


/  ,  (bk  .  \ 

-2kL  sin20Q  cos2  cos  0q 

a/T"-  sln20Q  sin2  <J>o  ^ 


(38) 


In  summary,  then,  for  incremental  strips  chosen  along  the  direction  of  propagation  In  the  half¬ 
plane  of  the  PO  current,  the  factor  (38)  should  multiply  the  fields  In  Eq.  (18)  when  Eq.  (18)  is  applied  to 
the  trailing  edge  of  the  truncated  half-plane. 

Of  course,  a  more  accurate  method  of  accounting  for  the  diffraction  at  the  tr  ailing  edge  of  the 
nonuniform  TE  current  emanating  from  the  leading  edge  of  an  arbitrary  scatterer,  would  be  to  first 
determine  the  path  of  the  diffracted  ray  from  the  leading  edge  to  the  trailing  edge,  and  then  calculate 
the  diffraction  at  the  trailing  edge  of  this  leading-edge  rav.  We  did  not  use  this  more  accurate 
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technique  in  the  present  report,  not  only  because  it  would  require  an  extra  calculation,  but  because  the 
TE  nonuniform  current  decays  fairly  rapidly  away  from  the  leading  edge  except  near  grazing  where  the 
diffracted  ray  path  nearly  coincides  with  the  direction  of  propagation  of  the  PO  current  in  the 
corresponding  half-plane  (and  thus  (38)  applies]. 

In  general,  Appendix  C  shows  that  whenever  the  significant  nonuniform  current  of  integration 
for  a  given  observation  point  extends  a  distance  from  the  edge  larger  than  a  small  fraction  (say  1/4)  of 
the  distance  to  the  focal  point  of  the  diffracted  rays,  the  incremental  strips  of  nonuniform  current  will 
not  closely  approximate  the  actual  nonuniform  current,  and  thus  the  fields  computed  from  the 
incremental  diffraction  coefficients  at  the  given  observation  point  may  contain  significant  error. 

(For  example,  see  the  discussion  of  Figures  22  -  25  in  Section  4.4.)  Near  grazing,  the  focal  length 
approaches  infinity  and  the  incremental  strips  of  nonuniform  current  accurately  approximate  the 
actual  nonuniform  current  across  the  whole  distance  from  leading  to  trailing  edges. 


2.6  Calculations  of  the  Two-Dimensional  Bistatic  Radar  Cross  Section  of  a  Strip 


In  this  section  we  show  the  results  of  calculations  made  of  the  two-dimensional  bistatic  cross 
section  of  an  infinitely  long  strip.  The  bistatic  cross  section,  o(<J>),  per  unit  length,  for  normal 
incidence  (that  is,  the  direction  vector  of  the  incident  plane  wave  normal  to  the  edges  of  the  strip)  of  an 
Infinitely  long  strip  is  defined  by  Eq.  (1.34)  of  Reference  12: 


®(4>)  =  p1J51«2;tp 


(39) 


where  Vs1  =  E®'1  if  the  electric  field  is  parallel  to  the  z-axis  (TM)  and  Vs-1  =  Hf-1  If  the  magnetic  field 
is  parallel  to  the  z-axis  (TE),  with  the  superscripts  s,i  denoting  the  scattered  and  Incident  field, 
respectively.  In  our  results  we  normalize  a  by  dividing  by  the  wavelength  X. 

In  Figures  3a  and  b  we  compare  the  TM  and  TE  back  scatter  cross  sections  obtained  using  the 
approximations  of  Sections  2.1,  2.3,  and  2.5,  with  the  corresponding  bistatic  cross  sections  obtained 
using  a  computer  program  to  calculate  the  exact  two-dimensional  fields  scattered  from  an  infinite 
strip. 13  The  width  of  the  strip  Is  s/X  =  30/2n:  =  4.77.  In  these  plots  0  is  the  angle  between  the  direction 
vector  of  the  incident  plane  wave  (or  the  observation  direction)  and  the  normal  to  the  plane  of  the 
strip.  For  calculating  the  contribution  of  the  leading  edge  of  the  strip  from  Eq.  (7a)  or  Eq.  (14a), 

<j>0  =  <J)  =  n/2  +  0,  while  for  calculating  the  contribution  of  the  trailing  edge  from  Eq.  (7a)  or  Eq.  (14a) 
with  Eq,  (35).  <J>0  -  4>  =  n/2  -  0,  0  <  0  <  n/2.  The  cross  sections  are,  of  course,  symmetrical  with  respect  to 
0  =  n/2.  The  total  field  scattered  from  the  strip  is  obtained  by  summing  the  contributions  of  the 
leading  and  trailing  edges.  Note  that  the  direction  of  z  In  Eqs.  (7a)  and  (14a)  must  be  reversed  for  the 
trailing  edge  from  what  it  is  for  the  leading  edge  of  the  strip.  As  can  be  seen  from  Figures  3a  and  3b,  the 

exact  back  scatter  cross  sections  ( - )  are  virtually  indistinguishable  from  the  approximate  back 

scatter  cross  sections  ( - ). 

In  Figures  4a  and  b  we  compare  the  approximate  TM  and  TE  specular  (or,  equivalently,  forward) 


13.  Dominek,  A.K.  (1988)  Personal  Communication,  Ohio  State  University  ElectroScience 
Laboratory. 
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Figure  3a.  Back  Scatt 


30. 


Figure  4b.  Specular  Scatter  Cross  Section  Pattern  of  Strip.  1 
- :  Exact. - :  PO  +  Approximate  Nonuniform  Ci 


scatter  cross  sections  with  the  corresponding  exact  cross  sections.  As  with  the  back  scatter  plots,  here 
0  is  the  angle  between  the  direction  vector  of  the  incident  plane  wave  (or  the  observation  direction)  and 
the  normal  to  the  plane  of  the  strip,  so  that  4>0  =  n/2  +  0  and  <|>  =  n/2  -  0  for  obtaining  the  contribution  of 
the  leading  edge  of  the  strip,  while  (j>0  =  n/2  -  0  and  <|>  =  n/2  +  0  for  the  trailing  edge.  Again,  there  is 
hardly  any  difference  between  the  exact  and  the  approximate  curves. 

In  Figures  5a  and  b,  we  compare  the  approximate  TM  and  TE  "side  scatter"  cross  sections  with  the 
corresponding  exact  cross  sections.  By  "side  scatter",  we  mean  that  the  observation  direction  and  the 
incidence  direction  are  at  right  angles  to  each  other.  In  these  plots  0,  is  the  angle  between  the  direction 
vector  of  the  incident  plane  wave  and  the  normal  to  the  plane  of  the  strip,  so  that  4>0  =  n/2  +  0,  and 
<J>  =  w,  foi  calculating  the  contribution  of  the  strip  leading  edge  from  Eq.  {7aj  or  Eq.  (14a),  while 
()>0  =  n/2  -  0,  and  $  =  x  -  0,  for  calculating  the  contribution  of  the  strip  trailing  edge  from  Eq.  (7a)  or 
from  Eqs.  (14a)  and  (35).  There  is  again  very  close  agreement  between  the  exact  and  the  approximate 
side  scatter  cross  section  curves. 

In  Figures  6a,  b,  and  c  we  compare  respectively  the  TE  back  scatter,  specular  (forward)  scatter, 
and  side  scatter  cross  section  plots  obtained  without  the  use  of  the  trailing  edge  correction  factor 
Eq  (35),  with  the  corresponding  exact  cross  section  plots.  Comparing  these  figures  with  Figures  3b,  4b, 
and  5b,  respectively,  we  see  that  the  trailing  edge  correction  factor  plays  a  significant  role  in 
improving  the  accuracy  of  the  approximation  in  the  back  scatter  and  side  scatter  cases. 

It  is  of  considerable  interest  to  see  the  limitations  of  the  PO  approximation  for  the  strip  bistatic 
cross  sections.  In  Figures  7a  and  b,  8a  and  b,  and  9a  and  b,  we  compare  the  TM  and  TE  back  scatter, 

specular  scatter,  and  side  scatter  cross  section  curves  obtained  with  the  PO  approximation  ( - ), 

with  the  exact  curves  ( - ).  It  Is  apparent  that  the  PO  approximation  is,  in  general,  of  rather 

limited  value  in  giving  an  accurate  picture  of  the  scattering  cross  section  patterns  apart  from  the 
portion  of  the  patterns  in  the  vicinity  of  specular  scatter. 

Figures  10a  and  b  compare  the  E-  and  H-wave  back  scatter  cross  sections  of  the  strip  computed 
from  the  high-frequency  solution  of  Michaeli's  1984  Radio  Science  paper14,  with  the  exact  E-  and 
H-wave  cross  sections.  Agreement  between  the  exact  E-wave  solution  and  Michaeli's  E-wave  solution 
is  excellent.  However,  Figure  10b  shows  a  significant  discrepancy  between  the  H-wave  solutions. 
Comparison  of  Figure  10b  with  Figure  6a  reveals  that  this  discrepancy  is  caused  by  the  neglect  in 
Michaeli's  solution  of  the  diffraction  of  the  H-wave,  leading-edge  nonuniform  current  at  the  trailing 
edge.  That  Is,  Michaeli's  high-frequency  solution  does  not  include  the  trailing  edge  factor  in  Eq.  (35), 
and  thus,  shows  poor  agreement  with  the  exact  solution  of  the  strip  illuminated  by  an  H-polanzed 
plane  wave. 

Ufimtsev’s  first  order  H-wave  back  scatter  radar  cross  section  from  his  1969  paper15  is  compareu 
to  the  exact  solution  in  Figure  11.  Except  in  the  vicinity  of  the  main  lobe,  this  first  order  H-wave 
solution  which  Is  more  accurate  (according  to  Uflmlsev15)  than  his  solution  in  Reference  1 1  used  by 
Coleman  et  al8,  shows  poor  agreement  with  the  exact  H-wave  solution. 


14  Michael!,  A  (1984)  A  closed  form  physical  theory  of  diffraction  solution  foi  electromagnetic 
scattering  by  strips  and  90°  dihedrals,  Radio  Science  19:609-616. 

15.  Ufimlsev,  P.  Ya.  (1969)  Asymptotic  investigation  of  the  problem  of  diffraction  on  a  strip.  Radio 
Engineering  and  Electronic  Physics  14:1014-1025. 
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Figure  5a.  Side  Scatter 
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Figure  6a.  Back 


26 


Figure  6b.  Specular  Scatter  Cross  Section  Pattern  of  Strip,  ks  =  30,  with  TE  (Parallel  Polarized)  Plane- Wave  Illumination; 
- :  Exact. - :  PO  +  Approximate  Nonuniform  Current  Field  (No  Trailing  Edge  Correction) 
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Figure  6c.  Side  Scatter  Cross  Section  Pattern  of  Strip,  ks  =  30.  with  TE  (Parallel  Polarized)  Plane-Wave  Illumination; 
- :  Exact. - :  PO  +  Approximate  Nonunifonn  Current  Field  (No  Trailing  Edge  Correction) 
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Figure  7a.  Back  Scatter  Cross  Section  Pattern  of  Strip,  ks  =  30,  with  TM  (Perpendicular  Polarized)  Plane-Wave  Illumination: 
- :  Exact, - :  PO 


Figure  8a.  Specular  Scatter  Cross  Section  Pattern  of  Strip,  ks  =  30,  with  TM  (Perpendicular  Polarized)  Plane-Wave 
Illumination; - :  Exact, - :  PO 
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Figure  10a.  Back  Scatter  Cross  Section  Pattern  of  Strip,  ks  =  SO,  with  TM  (Perpendicular  Polarized)  Plane-Wave  Illumination; 
- :  Exact, - :  Michaeli  (1984) 
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Scatter  Cross  Section  Pattern  of  Strip,  ks  —  30,  with  TE  (Parallel  Polarized)  Plane-Wave  Illumination; 
- ;  Ufimtsev  ( 1969) 


To  show  how  the  accuracy  of  the  approximations  derived  in  Sections  2.1  through  2.5  depends  on 
the  strip  width,  in  Figures  12a  and  b.  13a  and  b,  and  14a  and  b.  we  show  the  TM  and  TE  back  scatter, 
specular  scatter,  and  side  scatter  cross  section  plots  obtained  for  a  strip  of  width  s/X  =  3/2n  =  0.48. 
Although  the  approximations  are  not  quite  as  accurate  as  they  are  for  the  strip  of  width  s/X  =  30/2tc. 
nevertheless  the  approximate  cross  sections  differ  at  most  by  only  1.5  dB  from  the  exact  cross 
sections. 


3.  NONUNIFORM  (PTD)  INCREMENTAL  DIFFRACTION  COEFFICIENTS  FOR  THE 
TRUNCATED  HALF-PLANE 

In  this  section  we  obtain  the  TM  and  TE  nonuniform  or  PTD  IDC's  for  the  the  leading  and 
trailing  edges  of  a  truncated  half-plane.  The  TM  nonuniform  IDC  for  either  the  leading  or  trailing 
edge  of  a  truncated  half-plane  is  obtained  by  substituting  En0u^™'  in  Eq.  (24b)  of  Reference  3  with 
EnuCTM)  glven  by  coefficient  of  0g  in  Eq.  (10a).  Thus, 


jpnu(TM)/-  \  r  -><»  n  _ _  cin  _ 

dE  (r,)  ru,  -dz,Elz{  4jzr[  2V2  sin  2  sJn2 


<foi  sin  8t _ 1 


0OI  cos  4>ol  +  cos  a 


.  i  e~ln/4 


5(1-C°?i72  K°(ksSine 

1 1  -  cos  a  | 1/2  v 


0{  1 1  -  cos  a  I ) 


-2  cos  e-'ks  sln  e<» (cos  +  ccs  a)  F(2ks  sin  0O,  cos2 


V2"  cos^jfl-  e-,kss,n6oiwos$o{-  -os  aljU 


where 

sin  0j  cos  4>  j  +  cot  (cos  0(  +  cos  0oj)  ,  . 

cos  a  = - -7— a -  (41) 

sin  80j 

with  V|  the  angle  between  the  axis  of  the  incremental  strip  and  the  positive  local  z-axis.  We  have 
inserted  the  subscript  T  in  Eqs.  (40)  and  (41)  to  indicate  that  the  subscripted  quantities  are  defined 
with  reference  to  a  local  coordinate  system  with  origin  at  the  incremental  length  dz('. 

The  TE  nonuniform  IDC  for  the  leading  edge  of  a  truncated  half-plane  is  obtained  by  substituting 
H^nufTE)  ancj  pj^nufre)  ^  ^  (17)  0f  Reference  3  followed  by  Eq.  (18)  of  Reference  3,  with  Henu^  and 
H^nufTE)  given  hy  Eq.  (18a)  and  Eq.  (31a)  respectively.  Thus, 
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Specular  Scatter  Cross  Section  Pattern  of  Strip,  ks  =  3.  with  TM  (Perpendicular  Polarized)  Plane- Wave  Illumination; 
tact, - :  PO  +  Approximate  Nonuniform  Current  Field 
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Figure  13b.  Specular  Scatter  Cross  Section  Pattern  of  Strip,  ks  =  3,  with  TE  (Parallel  Polarized)  Plane-Wave  Illumination; 
_ _ — :  Exact, - :  PO  +  Approximate  Nonuniform  Current  Field 
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-60. 


o 


:  PO  +  Approximate  Nonuniform  Current  Field 


dEnurrc)(r,)  rUS~  -dzj  sign  (n  -  <J>0i)Z0H 


iz,  4itr{  sin  90( 


[  (  V2  e-'1^4  jVlkssln  0qi  (<»s«oj+co3a)  F^2ks  sin  0OI  cos2^j 


^2“l —  ^01 


V2  COS  —£• 

jT-'  COS  a  I  »/2  F°  ( ks  sin  eol  1 1  -  cos  a  I )  J 


+  [  1  -e 


-Iks  sin  Oqj  (cos  <>oi 


+  cos  a)  1  ^  cos  <(){  cos  9t  +  cos  a  cot  90t  sin  9; 
J  )  cos  <j)ot  +  cos  a 


(42) 


-f^2e-u/4  c~ikssineol(cos^l  +  cosa)  p  ^2ks  sin  9oI  cos2 


[*-" 


-1  ks  sin  0O|  (cos  <j>01  +  oos 


^COt  0O(  Sin  9t  J  9j 


-  cosC/cosa  e-'*/4  ('■“* "" <ta  h“4“*  “  F^ks  sin  e0,  cos*%l) 

V2  lcos|M  ,  r  1lA 

~]  i  -  cos  a |  »/2  F°(ks  sin 90ll  1  -cos  al)  J+  [  1  -e'lks  sln  Bab* hi* 006  “)  J  I 

with  cos  a  given  by  Eq.  (41).  The  TE  nonuniform  IDC  for  the  trailing  edge  of  a  trurK.dleu  naif-plane  is 
obtained  by  multiplying  the  edge  expression  in  Eq.  (42)  by  the  factor  given  in  Eq.  (35)  or  (37),  as 
explained  previously  in  Section  2.5. 


4.  SCATTERING  OF  A  PLANE  WAVE  BY  A  PERFECTLY  CONDUCTING  DISK 

In  this  section  we  apply  the  IDCs  derived  in  Section  3  to  calculate  the  far  fields  scattered  from  a 
perfectly  electrically  conducting  disk  of  radius  a.  The  geometry  is  shown  in  Figure  15.  The  primaiy 
source  is  a  plane  wave  whose  direction  of  propagation  lies  in  the  xz- plane  and  makes  an  angle  of  0, 
with  the  z-axis.  The  observation  direction  is  given  in  spherical  polar  coordinates  (R,  0,  <]>).  Two 
polarizations  of  the  incident  plane  wave  are  considered: 

1)  perpendicular  (TE)  or 

Finc=  Eie1K‘ ’  ry  .  (43a) 

Hinc  =  Y0E,e,k| '  r  (cos  0,  X  -  sin  0,  z  ,  (43b) 
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Figure  15.  Geometry  of  Disk  Illuminated  by  a  Plane  Wave 
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and 

2)  parallel  (TM)  or 

E |nc  =  E(  elk‘ '  r  (cos  9t  x  -  sin  9j  z)  .  (44a) 

Hinc=-Y0E1eiki  r'y  .  (44b) 

with 

k .  -  -kr ,  =  -k(sln  9(  x  +  cos  9j  z)  ,  (45) 

and 

—  A  A  A 

r  =  xx  +  yy  +  zz  .  (46) 

The  far  field  scattered  from  the  disk  is  obtained  by  adding  the  integrals  of  the  strip  nonuniform  IDC's 
around  the  edge  of  the  disk  to  the  PO  far  field.  The  nonuniform  high-frequency  solution  presented 
here  for  scattering  from  the  disk  remains  valid  for  all  angles  of  incidence  and  scattering,  and  thus,  is 
distinguished  from  the  previous  high-frequency  solutions  that  concentrated  on  back  scattering. 16'1'7 


4.1  Physical  Optics  Field 

The  PO  scattered  field  is  derived  in  Reference  18  and  we  give  only  the  results  here.  For 
perpendicular  polarization  of  the  incident  plane  wave 

ESI>0(P)  =  iE(a  ■■  —  p~— ■■  cos  9,(cos  9  sir.  $9  +  cos  <)><!>) 

while  for  parallel  polarization 

—  i>n  elkR  J i(kaP)  a  a 

Eg1  °(P)  =  tE,a  — ^ - p cos  9j(cos  9  cos  (j) 9  -  sin  $<{>) 

with 

(1  =  j^sin2  9  sin2  <}>  +  (sin  9S  +  sin  9  cos  <J>)2j'^2 
and  J  j  the  Bessel  function  of  order  one. 


(47) 


(48) 


(49) 


16  DeVore,  R.,  Hodge,  D.B.,  and  Kouyoumjlan,  R.C.  (1971)  Backscallering  cross  sections  of  circular 

disks  for  arbitraiy  incidence,  J.  Applied  Physics  42:3075-3083. 

17  Maryland,  D.P.,  Balanis,  C.A.,  and  Brumley,  3.A.  (1987)  Higher  order  diffractions  irom  a 

circular  disk,  IEEE  Trans.  Antennas  Propagat.  AP-35: 1436-1444. 

13.  Troll,  K.  (1988)  The  Disk:  A  Comparison  of  Electromagnetic  Scattering  Solutions  and  Its  Use  as 
a  Calibration  Standard J'oi  Bistatic  RCS  Measurements,  RADC-TR-88-16,  ADA200327. 


46 


4.2  Transformation  of  the  Local  Coordinates  of  the  Strip  IDC's  to  the  Global 
Coordinate  System  of  the  Disk 


The  correction  to  the  PO  electric  far  field  of  the  disk  is  obtained  as  indicated  above  by  integrating 
the  nonuniform  strip  IDC's.  Eqs.  (40)  and  (42),  around  the  edge  of  the  disk.  Since  all  quantities  in 
Eqs.  (40)  and  (42)  are  defined  with  respect  to  a  local  coordinate  system  with  origin  at  the  differential 
element  of  the  disk  edge,  it  is  first  necessary  to  transform  Eqs.  (40)  and  (42)  to  the  global  coordinate 
system  of  Figure  15.  It  is  also  necessary  to  define  the  local  coordinate  system.  Since  we  are  modelling 
the  edge  locally  by  the  truncated  half- plane  defined  as  in  Figure  2,  we  want  the  local  x-axls  to  be  in  the 
plane  of  the  disk,  normal  to  the  edge,  and  directed  inward  from  the  edge;  the  local  y-axis  to  be  the 
normal  to  the  disk  at  the  edge  directed  toward  the  half-space  from  which  the  illuminating  plane  wave 
is  incident;  and  the  local  z-axis  to  be  tangent  to  the  edge  of  the  disk.  Hence, 


A  A  |  i  A  i  A 

x,  = -p  =  -  cos  4>  x  -  sin  «J>  y  , 


(50a) 


A  A 

y i~  7-  < 


(50b) 


and 


A  A  A  A  ,  ,  A  ,  A 

z I  =  x j  x  y  (  =  q  =  cos  <J>  x  -  sin  4>  y  . 


(50c) 


where  we  use  primes  to  denote  the  point  of  integration. 

We  now  systematically  express  the  locally  defined  quantities  of  Eqs.  (40)  and  (42)  in  the  global 
coordinate  system.  Starting  with  0OI,  the  angle  between  the  local  unit  z-vector,  4>\  and  the  vector  r, 
given  by  Eq.  (45), 

A  A 

cos  0O,  =  r ,  •  z  ( 


=  -  sin  0,  sin  4>'  (51a) 

and 

sin  0ol  =  +  (1  -  sin2  0,  sin2  $') 1/2  (51b) 

where  the  positive  root  must  be  taken  since  0  <  0O!  £  x. 

Next,  since  4>0I  is  the  angle  between  x(  and  the  projection  of  r,  on  the  X|  y(  -plane, 

cos  o|  =  r,  -  x(  /sin  0o, 

=  -sin  0j  cos  <J> ' /sin  0O{  (52a) 

and 

sin  4> o|  =ri-yl  /sin  0ol 

=  cos  0( /sin  0O(  (52b) 

from  which 
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cos  2~  =  + 


sinGjCOS^' M/2 
sin  0oj 

2  )  • 


where  the  positive  root  must  be  taken  since  for  a  plane  wave  incident  from  the  upper  half-space,  z  >  0. 
0<<t>ol<7t. 

Proceeding  to  0(,  the  angle  between  z(  and  the  ray  from  the  edge  point  to  the  far-fleld  observation 
point, 

A  A  A,  A 

cos  0{  =  Z|  •  r  (  =  <}>  •  R 


=  -sin  0  sin(<S>'  -  4>) 


sin  0j  =  +  £  1  -  sin2  0  sin2 (<)>'  -  <J> )  j^2 

where  the  positive  root  must  be  taken  since  0  <  0[  £  n.  The  cosine  and  sine  of  are  given  by 

A  A  A  A 

cos  <(>i  =  r  |  •  xj  /sin  0j  =  R  •  x{  /sin  0( 

=  -sin  0cos(4>'  -  4> ) /sin  0t 


A  A  ,  A  A  , 

sin  (J)  |  =  r  (  •  y ;  /sin  0t  =  R  •  y  t  /sin  0; 

=  cos  0/sin  0,  . 

The  unit  vector  0  ( is  transformed  to  global  coordinates  by  starting  with 

A  A  A  A 

0|  -■  cos  Oj  cos  i)>(  Xj  +  cos  0/  sin  (j)t  y  j  -  sin  0j  z(  , 

substituting  global  Cartesian  component  expressions  for  Xj,  yt ,  and  z(  ,  and  simplifying 
trigonometrically,  thereby  obtaining 

0 1  =  (-cos  0j  cos  4>(  cos  $'  +  sin  0;  sin  <J>')x 
-  (cos  Oj  cos  j  sin  ij) '  +  sin  0;  cos  4> '  )y 


+  cos  0;  sin  <j)(  z  ,  (56) 

with  cos  0( ,  sin  0j ,  cos  <}>; ,  and  sin  (j), ,  given  by  Eqs.  (54a),  (54b),  (55a),  and  (55b)  respectively.  A  similar 
procedure  gives 

f  -  sin  <{>  j  cos  ij>  x  +  sin  <p  j  sin  <p  y  +  cos  (pj  z  .  (57) 

When  the  IDC's,  Eqs.  (40)  and  (42),  are  integrated  around  the  edge  of  the  disk  to  obtain  the 
nonuniform  current  correction  to  the  PO  far  field,  it  is  the  0-  and  <f>-  components  of  the  electric  far  field 
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that  are  required  (the  radial  component,  is  of  course,  zero).  Hence,  it  is  desirable  to  have  expressions 


A  A 

for  the  0-  and  <J>-  components  of  0  ( and  <J> ; .  Letting 

AA  AA  AA  AA  . 

cr  =  0  j  •  0  ,  r  =  0;-<t>,  p  =  4>  i  •  0  ,  v  =  <)>  {  •  <(>  ,  (58) 

we  have 

a  =  a  cos  0  cos  <J>  +  b  cos  0  sin  <j>  -  c  sin  0  ,  (59a) 

x  =  -a  sin  <(>  +  b  cos  <(»,  (59b) 

p  =  d  cos  0  cos  <}>  +  e  cos  0  sin  <j>  -  f  sin  0  ,  (59c) 


V  =  -d  sin  y  i  e  cos  6.  (59H) 

where  a,  b,  and  c  are  the  coefficients  of  x,  y ,  and  z .  respectively,  in  Eq.  (56)  and  d,  e,  and  f  are  the 
coefficients  of  x ,  y ,  and  z  in  Eq.  (57).  Substituting  explicit  expressions  for  a,  b,  and  c,  and  d,  e,  and  f 
from  Eqs  (56)  and  (57)  in  (59)  followed  by  substitution  of  the  expressions  Eqs.  (54)  and  (55)  for  cos  0( . 
sin  0, ,  cos  4>t .  and  sin  <J>; .  yields  after  algebraic  and  trigonometric  simplification 


_ cos  0  sin  u _ 

°  _  (1  -  sin2  0  sin2  u)  '/2 


(60a) 


_ cos  u _ 

x  “  "(1  -  sin2  0  sin2  u)1/2 


(60b) 


_ cos  u _ 

“  (1  -  sin2  0  sin2  u)1'2  ’ 

_ cos  0  sin  u 

V  ""  ( 1  -  sin2  0  sin2  u) 1//2 

where 

u  =  <J>'  -  <|>  . 

Next,  proceeding  to  elkri/ri  and  replacing  rtby  r 
eikr;  elkr  e i k ( R  —  r’- R) 

r, . =  r  R  ■ 

It  is  then  simple  to  find  that 

r  ’  •  R  =  a  sin  0  cos(4> '  -  40  • 

The  differential  length  dz’,  Is  given  by 
dz'j  =  ad* ' 


(60c) 

(60d) 


(61) 


(62) 


(63) 


(64) 


Finally,  the  components  of  the  illuminating  field  in  the  local  z-directlon,  E ^  and  Hlzi,  are  given  by 
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(65a) 


EiZ|  -  Ei«'  • 


H„,=  I V*  (65b) 

where,  from  Eqs.  (43)  and  (44), 

E,^  =  Etcos<t)'  elEi-f‘  ,  (66a) 

=  -Y0EjCOS  0j  sin  <J» '  e*  *9 ' r '  ,  (66b) 

for  perpendicular  polarization  of  the  incident  plane  wave,  and 

£(<■,■  =  -E,cos  0t  sin  <J) '  e*kt ' r  ' ,  (67a) 

H^^-YoE.cos-C'  (67b) 

for  parallel  polarization  of  the  incident  plane  wave. 

4.3  Choice  of  Incremental  Strips 

In  integrating  the  nonuniform  current  IDC's  around  the  edge  of  the  disk  to  obtain  the  correction 
to  the  PO  far  field,  a  choice  must  be  made  of  the  orientation  of  the  incremental  strips,  since  the  strips 
can  be  chosen  to  cut  across  the  disk  in  any  direction.  Thus,  for  example,  the  strips  can  be  taken 
straight  across  the  disk  parallel  to  the  x-axls,  the  projection  of  the  direction  vector  of  the  incident 
plane  wave  on  the  disk,  or  they  can  be  chosen  to  lie  in  the  directions  of  the  diffracted  rays  on  the  disk, 
or  they  can  be  chosen  to  be  diameters,  normal  to  the  edge  at  the  points  of  integration.  In  general,  if  the 
nonuniform  current  dies  away  rapidly  from  the  edge,  the  choice  of  the  orientation  of  the  incremental 
strips  plays  little  role  in  determining  the  nonuniform  current  field.  However,  if  the  nonuniform 
current  does  not  decrease  rapidly  away  from  the  edge  as  it  does  not  here,  for  parallel  polarization  close 
to  grazing  incidence,  the  choice  of  the  incremental  strips  can  significantly  affect  the  calculated 
nonuniform  current  field. 

In  the  calculations  performed,  two  choices  of  incremental  strips  were  made;  1)  straight  across 
the  disk  parallel  to  the  x-axis;  and  2)  In  the  directions  of  the  diffracted  rays  on  the  disk.  (The  choice  of 
diameter  incremental  strips  was  ruled  out  because  it  would  lead  to  serious  overlapping,  and  hence 
overweighting,  of  the  nonuniform  currents  close  to  the  center  of  the  disk.)  The  choice  of  the 
incremental  strip  orientations  enters  into  the  calculations  analytically  through  1)  the  value  of  s,  the 
width  of  the  strip,  in  the  IDC's  given  by  Eqs.  (40)  and  (42);  and  2)  cot  \j/j  in  Eq.  (41)  where  is  the  angle 
between  the  axis  of  the  incremental  strip  and  the  positive  local  z-axis.  The  derivation  of  these 
quantities  is  given  in  Appendix  B  and  we  state  the  results  here.  For  incremental  strips  taken  straight 
across  the  disk  parallel  to  the  x-axis 


s  =  2a  cos2  (J)’  , 

(68) 

cot  \j/(  =  tan  <j> '  , 

(69) 
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(70) 


while  for  strips  taken  along  the  diffracted  rays  on  the  disk 
s  =  2a(  1  -  sin2  0|  sin2  <)>'), 


cot  =  -cos  0oi  /sin  0O(  .  (71) 

with  cos  0O1  and  sin  0O1  given  by  Eq.  (51).  For  both  choices,  the  parameter  L  in  the  factors,  (37)  or  (38), 
used  to  multiply  the  TE  trailing  edge  field,  is  taken  as  the  value  for  strips  straight  across  the  disk;  that 
is. 

L  —  2a  !  CuS  4-  i  (72) 


because  the  use  of  (37)  as  a  simple  multiplicative  factor  applied  to  the  incident  field  assumes  that  the 
direction  of  propagation  of  the  field  impinging  upon  the  trailing  edge  is  In  the  direction  of  the  incident 
field.  Near  grazing,  where  the  multiplicative  factor,  (37),  differs  significantly  from  unity,  all  the 
diffracted  rays  approach  the  direction  straight  across  the  disk  and  thus,  the  value  of  L  in  Eq.  (72)  is 
valid  for  the  diffracted  ray  directions  as  well  as  the  direction  across  the  disk.  Away  from  grazing, 
where  the  diffracted  rays  do  not  all  lie  in  the  direction  across  the  disk,  the  multiplicative  factor 
approaches  unity,  and  thus  its  effect  becomes  insignificant  for  large  disks. 

4.4  Calculations  of  the  Bistatic  Radar  Cross  Section  of  a  Disk 

In  this  section  we  show  the  results  of  calculations  of  the  bistatic  section  of  a  perfectly  conducting 
disk.  The  bistatic  cross  section.  c(0,<J>),  is  defined  by  Eq.  (1.30)  of  Reference  12: 


a(0,4>) 


11  m 
R-  >oo 


4xR2 


IesI2 


(73) 


where  E  s  is  the  scattered  electric  field  at  the  observation  point  (R,0,<]>).  In  our  calculations  we 
normalize  the  bistatic  cross  section  by  dividing  it  by  the  square  of  the  wavelength.  We  present  plots  of 
the  blstatic  cross  section  for  perpendicular  and  parallel  polarization  of  the  Incident  plane  wave,  and 
for  both  co-  and  cross-polarized  scattered  fields  defined  according  to  Ludwig's  third  definition19  as 


=  E0cos  4>  +  E<j,  sin  4>  , 

(74a) 

=  E0cos  <}>  —  E^  sin  4>  - 

(74b) 

Cross  sections  were  calculated  for  the  two  choices  of  Incremental  strip  orientations  discussed  in  the 
previous  section;  strips  taken  straight  across  the  disk  parallel  to  the  x-axis,  the  projection  of  the 
direction  vector  of  the  incident  plane  wave  on  the  disk;  and  strips  taken  in  the  direction  of  the 
diffracted  rays  on  the  disk.  The  approximate  cross  sections  obtained  by  adding  the  integrals  of  the 
nonuniform  current  IDC's  around  the  edge  of  the  disk  to  the  PO  fields  are  compared  with  the  exact 


19.  Ludwig,  A.C.  (1973)  The  definition  of  cross  polarization.  IEEE  Trans.  Antennas  FTopagat. 
AP-21:116-I19. 
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cross  sections  obtained  from  an  eigenfunction  solution.20  A  modification  of  a  computer  code 
implementation20  of  the  eigenfunction  solution,  made  by  Dominek21  to  allow  for  the  use  of  larger  ka 
values,  was  used  to  calculate  the  exact  cross  sections. 

In  Figures  16a  and  b,  we  show  the  back  scatter  (0  =  0,,  <}>  =  <J>,  =  0)  cross  section  patterns  for  a  disk  of 
size  ka  =  15  for  perpendicular  polarization  of  the  incident  plane  wave,  with  the  incremental  strips 
taken  straight  across  the  disk  and  in  the  direction  of  the  diffracted  rays,  respectively.  Figures  17a  and 
b  show  the  corresponding  back  scatter  patterns  for  parallel  polarization  of  the  incident  plane  wave. 
Both  strip  orientations  yield  similar  satisfactory  approximations  over  the  entire  range  of  0. 

In  Figures  18a  and  b  we  show  the  specular  scatter  (0  =  0,,  <J>  =  n)  cross  sections  for  perpendicular 
polarization  of  the  incident  plane  wave,  with  the  strips  oriented  parallel  to  the  x-ax»«  and  in  the 
direction  uf  the  diffrueted  iu)a,  leopcelivcly.  Figures  19a  and  b  show  the  corresponding  plots  for 
parallel  incident  plane  wave  polarization.  Both  choices  of  incremental  strips  yield  excellent 
approximations  over  the  entire  range  of  0  that  could  be  used  as  an  accurate  calibration  standard  for 
bistatic  RCS  measurements. 

The  side  scatter  (observation  direction  normal  to  the  incidence  direction:  0  =  n/2  —  01,  =  n)  cross 

section  patterns  for  perpendicular  polarization  of  the  incident  plane  wave,  with  the  incremental 
strips  parallel  to  the  x  axis  and  in  the  directions  of  the  diffracted  rays,  are  shown  in  Figures  20a  and  b, 
respectively.  The  corresponding  patterns  for  parallel  incident  plane  wave  polarization  are  shown  in 
Figures  21a  and  b.  The  patterns  obtained  with  the  strips  taken  parallel  to  the  x-axis  give  a  slightly 
better  approximation  to  the  exact  cross  sections  than  do  the  patterns  obtained  with  the  strips  in  the 
direction  of  the  diffracted  rays. 

A  significant  feature  seen  in  Figures  17,  19,  and  21  is  that  the  scattered  fields  of  the  high- 
frequency  solution  approach  zero  like  the  exact  scattered  fields  as  the  parallel  polarized  plane  wave 
approaches  grazing  incidence.  This  high  accuracy  near  grazing  for  the  parallel  polarized  plane  wave 
further  confirms  the  validity  of  the  trailing  edge  factor  in  (37). 

In  Figures  22  -  25  we  show  the  cut  through  the  cross  section  pattern  in  the  plane  defined  by  4)  =  45° 
and  <{>  =  225°,  for  a  plane  wave  incident  on  the  disk  at  angle  of  0,  =  45°.  In  these  plots  we  have  let 
Jiy  =  90°-  0  for  <J>  =  45°,  90°  >  0  >  0°.  and  y  =  90°  +  0  for  <j>  =  225°,  0  <  0  <  90°.  We  show  both  the  co-  and  cross¬ 
polarization  patterns  in  this  cut.  (For  the  back,  specular,  and  side  scattei  patterns,  the  cross-polarized 
field  equals  zero.)  Figures  22a  and  b  display  the  co- polarized  pattern  for  perpendicular  polarization  of 
the  incident  plane  wave,  with  the  incremental  strips  parallel  to  the  x-axis,  and  in  the  diffracted  ray 
directions  respectively.  Figures  23a  and  b  show  the  corresponding  cross-polarized  patterns.  In 
Figures  24a  and  b  we  show  the  co-polarized  patterns  for  parallel  polarization  of  the  incident  plane 
wave,  and  in  Figures  25a  and  b  the  corresponding  cross-polarized  patterns  are  shown. 

We  see  from  Figures  22-25  that  again  there  is  not  a  large  difference  in  accuracy  between  the 
high  frequency  solution  calculated  from  strips  parallel  to  the  x-axis  and  in  the  directions  of  the 
diffracted  rays.  The  larger  discrepancy  between  the  high-frequency  solution  and  the  exact  solution  for 
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Figure  16a.  Back  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Perpendicular  Polarized  Plane-Wave  Illumination  an< 

Incremental  Strips  Parallel  to  x-Axis;  - :  Exact, - :  PO  +  Approximate  Nonuniform  Current  Field  Obtained  by 

Integrating  IDC 
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Figure  16b.  Back  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Perpendicular  Polarized  Plane-Wave  Illumination  and 

Incremental  Strips  in  Direction  of  Diffracted  Rays:  - :  Exact. - :  PO  +  Approximate  Nonuniform  Current  Field 

Obtained  by  Integrating  IDC 


Figure  17a.  Back  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Parallel  Polarized  Plane-Wave  Illumination  and 

Incremental  Strips  Parallel  to  x-Axis:  - :  Exact, - :  PO  +  Approximate  Nonuniform  Current  Field  Obtained  by 

Integrating  IDC 
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Figure  I9a.  Specular  Scatter  Cross  Section  Pattern  cf  Disk,  ka  =  15.  with  Parallel  Polarized  Plane- Wave  Illumination  and 

.ncremental  Strips  Parallel  to  x-Axis;  - :  Exact, - :  PO  +  Approximate  Nonuniform  Current  Field  Obtained  by 

Integrating  IDC 
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Figure  21b.  Side  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Parallel  Polarized  Plane-Wave  Illumination  and 

Incremental  Strips  in  Direction  of  Diffracted  Rays:  - :  Exact. - :  PO  +  Approximate  Nonuniform  Current  Field 

Obtained  by  Integrating  IDC 
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Figure  22b.  Co-Polarized  Cross  Section  Pattern  of  Disk,  ka  =  15,  Illuminated  by  Perpendicular  Polarized  Plane  Wave  at 

0,  =  45°  in  Plane  Defined  by  <+>  =  45°  and  <■>  =  225°,  with  Incremental  Strips  in  Direction  of  Diffracted  Rays:  - : 

Exact, - :  PO  +  Approximate  Nonuniform  Current  Field  Obtained  by  Integrating  IDC 
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Figure  24a.  Co-Polarized  Cross  Section  Pattern  of  Disk,  ka  =  15,  Illuminated  by  Parallel  Polarized 

Plane  Defined  by  b  =  45°  and  6  =  225°,  with  Incremental  Strips  Parallel  to  x-Axis;  - :  Exact 

Approximate  Nonuniform  Current  Field  Obtained  by  Integrating  IDC 


e  24b.  Co-Polarized  Cross  Section  Pattern  of  Disk,  ka  =  15.  Illuminated  by  Parallel  Polarized  Plane  Wave  at  9,  =  45' 

;  Defined  by  d  =  45°  and  b  =  225°,  with  Incremental  Strips  in  Direction  of  Diffracted  Rays;  - :  Exact. - 

Approximate  Nonuniform  Current  Field  Obtained  by  Integrating  IDC 


these  45°  patterns  are  caused  by  significant  secondary  diffraction  between  the  leading  and  trailing 
edges  that  is  not  accurately  accounted  for  away  from  grazing  incident  by  the  incremental  strip 
approximation  for  the  current  on  this  disk  of  fairly  small  radius  a  =  15/2it  X  «  2AX.  This  discrepancy 
diminishes  for  electrically  larger  disks.  (See  Appendix  C.) 

As  with  the  strip,  it  is  of  interest  to  see  the  limitations  ci  i  ,  PO  approximation  for  the  disk 

bistatic  cross  sections.  In  Figures  26a  and  b,  27a  and  b,  and  28a  and  b,  we  compare  the  PO  ( - )  and 

exact  ( - )  back  scatter,  specular  scatter,  and  side  scatter  cross  section  patterns  for  perpendicular 

and  parallel  polarization  of  the  illuminating  plane  wave.  Figures  29a  and  b  compare  the  co-and 
cross-polarized  PO  and  exact  patterns  In  the  plane  defined  by  <}>  =  45°  and  y  =  225’  for  a  perpendicular 
polarized  plane  wave  incident  on  the  disk  at  an  angle  of  0f  =  45°.  Figures  30a  and  b  show  the 
corresponding  patterns  for  parallel  polarization  of  the  incident  plane  wave.  It  is  apparent  from  these 
figures  that  the  PO  approximation  does  not.  In  general,  give  an  accurate  representation  of  the  bistatic 
cross  section  patterns,  especially  in  the  vicinity  of  pattern  minima.  This,  of  course,  serves  to 
emphasize  the  importance  of  accounting  for  the  nonuniform  current  field  in  calculating  radar  cross 
sections  of  objects  with  edge  discontinuities. 

To  sec  how  the  accuracy  of  our  approximation  depends  on  the  disk  size,  in  Figures  31a  and  b,  32a 
and  b,  and  33a  and  b.  we  show  the  perpendicular  and  parallel  polarization  back  scatter,  specular 
scatter,  and  side  scatter  patterns  for  a  disk  of  size  ka  =  1.5.  For  these  patterns  the  incremental  strips 
have  been  taken  parallel  to  the  X -axis.  Although  there  are  differences  of  a  few  dB  between  the  exact 
and  IDC  solutions  over  much  of  the  range  of  0,  the  pattern  shapes  are  surprisingly  accurate  for  this 
u.mall  disk  of  diameter  less  than  one-half  a  wavelength. 
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Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Perpendicular  Polarized  Plane-Wave  Illumination; 
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Figure  27a.  Specular 
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Figure  27b.  Specular  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Parallel  Polarized  Plane-Wave  Illumination: 
- :  Exact, - :  PO 
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Figure  28a.  Side  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Perpendicular  Polarized  Plane-Wave  Illumination: 
- :  Exact, - :  PO 
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Figure  28b.  Side  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  15,  with  Parallel  Polarized  Plane-Wave  Illumina 
- :  Exact, - :  PO 
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Figure  Ola.  Back  Scatter  Cross  Section  Pattern  of  Disk,  ka  =  1.5,  with  Perpendicular  Polarized  Plane-Wave  Illumination  and 

Incremental  Strips  Parallel  to  x-Axis;  - :  Exact, - :  PO  +  Approximate  Nonuniform  Current  Field  Obtained  by 

Integrating  IDC 
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Appendix  A:  Derivation  of  an  Integral 
of  the  Fresnel  Integral 


In  this  Appendix  we  obtain  a  closed  form  expression  for  the  following  integral  of  the  Fresnel 
integral: 

X 

j  F(u)elaudu  .x  >  0  .  (A-l) 

o 

where  F(u)  is  the  Fresnel  integral  defined  by 

U 

1  C 

F(u)  =-==  ~=  dt  (A- 2) 

v2jc  J  Vt 
o 

To  integrate  (A-l),  begin  by  integrating  by  parts  to  obtain 

1  r  elta*1)u 

J  ^[u 

o 

We  now  consider  the  following  three  cases:  i)  a  +  1  >  0;  ii)  a  +  1  <  0:  and  ill)  a  +  1  =  0. 
i)  a  +  1  >  0 

For  a  +  1  >  0  let  t  =  (a  +  l)u.  Then  from  Eq.  (A-3) 
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11)  a+  1  <0 

Let  t  =  la  +  1  |u.  Then  from  Eq.  (A-3) 
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Hi)  a  +  1  =  0  or  a  =  -1 


From  Eq.  (A-3) 
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Appendix  B:  Derivation  of  Incremental 

Strip  Parameters 


In  this  Appendix  we  derive  expressions  for  the  width,  s,  of  the  incremental  strip  in  the  IDC’s 
given  by  Eqs.  (40)  and  (42),  and  for  cot  \|/j  in  Eq.  (41)  where  \y,  is  the  angle  between  the  axis  of  the 
incremental  strip  and  the  positive  local  z-direction.  Both  quantities  depend  on  the  choice  of 
orientation  of  the  incremental  strips.  We  consider  two  choices:  1)  incremental  strips  chosen  parallel 
to  the  x-axis  (the  projection  on  the  plane  of  the  disk  of  the  direction  vector  of  the  incident  plane  wave); 
and  2)  incremental  strips  chosen  in  the  direction  of  the  diffracted  rays  on  the  disk. 

Bl.  Incremental  Strips  Parallel  to  the  X-Axis 

Referring  to  Figure  Bla,  we  consider  an  incremental  strip  taken  parallel  to  the  x-axis  at  the  point 
on  the  edge  of  the  disk  at  4>',  -rc  <  4>'  <  n.  Since  the  local  unit  z-vector,  z(,  is  equal  to  [(Eq.  (50c)],  <j>j  is 
the  angle  between  the  chord  parallel  to  the  x-axis  at  4>'  and  the  tangent  to  the  disk  at  <J>'.  From 
elementary  geometric  considerations 
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so  that 

cot  vj/(  =  tan  4> *  . 


(B2) 


Figure  Bla.  Geomeliy  of  Disk  with 
Incremental  Strips  Parallel  to  the  x-Axis 


Figure  Bib.  Corresponding 
Strip  Geometry 


Since  the  disk  is  modelled  locally  by  a  strip  whose  edges  are  parallel  to  the  local  z-axis  (see  Figure 
Bib),  the  strip  width  s  In  the  IDC's  Eqs.  (40)  and  (42)  is  related  to  the  chord  length  L  in  Figure  Bla  by 

s  =  L  sin  .  ( 

Hence,  with  Eq.  (Bl) 

s  =  L  I  cos  <J) '  I  .  ( 


L  =  2a  I  cos  <}) '  I 


so  that 


s  =  2a  cos2  (t>‘ 


B2,  Incremental  Strips  in  the  Direction  of  the  Diffracted  Rays 


A  A  A 

rd  =  ax  +  Py 


a2  +  p2  =  1  (E 

be  (he  diffracted  rav  unit  vector  on  the  disk  at  the  point  (x'.  y')  =  (a  cos  <{>',  a  sin  <}>').  The  angle.  that  r 

A, 

makes  with  the  local  unit  z-vector,  4>  is  given  by 

l)/j  =  71  -  0ot  (E 


where  0oi  is  the  angle  between  the  local  z-axis  and  the  vector  in  the  direction  from  which  the 
illuminating  plane  wave  is  incident.  Then  from  Eqs.  (51a,b) 

cos  =  -  cos  G0(  =  sin  04  sin  <J> 1  ,  (BlOa) 

sin  \|/;  =  sin  0oi  =  (1  -  sin2  0{  sin2  t}> ') 1  /2  ,  (BlOb) 


so  that 

cot  Vj  =  -  cot  0oj  .  (Bll) 

To  determine  the  length  of  the  chord,  L,  made  by  the  diffracted  ray  on  the  disk  we  first  need 
expressions  for  the  direction  cosines  a  and  P  in  Eq.  (B7). 

Since 

rd  •  <}> '  =  cos  \|/(  {B 12) 

it  follows  from  Eq.  (BlOa)  that 

-a  sin  1  +  p  cos  $ 1  =  sin  0t  sin  (j>'  (B13) 

or 


P  =  tan  <{) ' (oc  +  sin  0t)  .  (B14) 

Substituting  Eq.  (B14)  into  (B8),  solving  the  resulting  quadratic  equation  for  a,  and  substituting  a  in 
Eq.  (B14)  gives 

a  =  -sin  0(sin2  <!>'  ±  cos  <t>'(l  -  sin2  0i  sin2  (jO1/2  ,  (B15a) 

P  =  sin  ({t'J^sin  0,  cos  <J)'  ±  (1  -  sin2  0j  sin2  <)> ' } 1  /2J  .  (B15b) 

Now  let 

x  =  a  cos  <t>'  +  au  ,  (B16a) 

y  =  a  sin  «j» '  +  Pu  ,  (B16b) 

be  the  parametric  representation  of  points  on  the  diffracted  ray.  Substituting  for  a  and  p  from  Eq. 
(B15)  inEq.  (B 16).  setting 

x2  +  y2  =  a2  (B17) 

and  letting  L  be  the  resulting  solution  for  u,  yields 

L  =  -2a  •  ±(  1  -  sin2 0, sin2 (J> ') 1  /2  (B18) 

where  the  ±  corresponds  to  that  in  Eq.  (B15).  It  follows  that  the  minus  sign  must  be  chosen  in 
Eq.  (B15)  and 

L  =  2a  (1  -  sin2  0,  sin2  $’)  ‘/2  (B19) 

Hcncc,  using  Eqs.  (B3)  and  (BiOb), 

s  =  2a  ( 1  -  sin2  0j  sin2  <J)')  .  (B20) 
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Appendix  C:  Conditions  for  Incremental  Current 
Strips  to  Accurately  Predict  Fields 
Diffracted  from  the  Trailing  Edge 


In  Section  2.5  we  modified  the  TE  fields  Incident  upon  the  trailing  edge  of  a  truncated  half-plane 
to  account  for  the  nonuniform  TE  current  that  emanates  from  the  leading  edge  and  diffracts  at  the 
trailing  edge.  (The  nonuniform  TM  currents  decay  rapidly  away  from  the  leading  edge  for  all  angles  of 
Incidence,  and  thus,  produce  negligible  diffracted  fields  at  the  trailing  edge  of  a  truncated  half-plane  a 
few  or  more  wavelengths  across.)  In  applying  the  TE  modification  derived  in  Section  2.5  to  the  disk, 
we  assume  that  the  incremental  nonuniform  current  strips  mail  .tain  a  constant  width  across  the 
disk.  Although  this  assumption  is  valid  for  the  infinitely  long  straight  edges  of  a  truncated  half¬ 
plane,  it  may  not  yield  an  accurate  representation  of  the  actual  nonuniform  currents  beyond  a  certain 
distance  from  the  leading  edge  of  a  finite  dimensional  flat  plate  such  as  the  disk.  The  purpose  of  this 
Appendix  is  to  establish  the  distance  from  the  leading  edge  of  a  flat  plate  that  the  nonuniform  current 
predicted  by  incremental  half- plane  current  strips  remains  a  good  approximation  to  the  actual 
current. 

When  the  incident  plane-wave  is  grazing,  the  magnitude  of  the  nonuniform  TE  current  remains 
constant  across  the  flat  plate  and  the  nonuniform  incremental-strip  current  agrees  exactly  with  the 
ac‘ual  current.  For  an  incident  plane-wave  far  from  grazing,  the  actual  TE  nonuniform  current, 
^actual*  emanating  from  the  leading  edge  of  the  plate  can  be  evaluated  approximately  by  a  stationary 
phase  evaluation  of  the  diffraction  integral.  This  evaluation  leads  to  the  usual  geometric  factor  In  the 
edge-diffracted  rayscl;  specifically 


Cl.  Kouyoumjian,  R.G.  (1975)  The  geometrical  theory  of  diffraction  and  its  application,  Chapter  6  in 
Vol.  3  of  Topics  in  Applied  Physics  (Numerical  and  Asymptotic  Techniques  in 
Electromagnetics),  R.  Millra,  Ed.,  New  York:  Springer  Verlag. 
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where  d  is  the  distance  from  the  leading  edge  (diffraction  point)  along  the  diffracted  ray  to  the 
observation  point  in  the  plate,  and  q  is  the  focal  length  along  this  diffracted  ray  measured  from  the 
leading  edge.  The  focal  length  is  determined  by  the  curvature  of  the  leading  edge  at  the  diffraction 
point  and  the  direction  of  incidence  of  the  Illuminating  plane  wave.  The  constant  A  is  a  complex 
amplitude  proportional  to  the  amplitude  of  the  incident  plane  wave.  The  derivation  of  Eq.  (Cl) 
assumes  that  the  distance  d  from  the  leading  edge  to  the  observation  point  is  greater  than  about  a 
wavelength,  and  that  the  incident  plane-wave  is  sufficiently  far  from  grazing  for  the  nonuniform  TE 
current  to  decay  as  1/Vd  for  \  <  d  «  q. 

The  nonuniform  TE  current,  predicted  by  incremental  half-plane  current  strips 

of  constant  width,  taken  along  the  directions  of  the  diffracted  rays  through  the  trailing  edge  points,  is 
given  approximately  by 

(C2) 

for  d  >  X  and  the  incident  plane- wave  far  from  grazing.  The  1  /  Vd  factor  can  be  derived  by  subtracting 
the  PO  current  from  the  total  current  in  Eq.  (32a)  and  then  letting  2k  sin  80  cos2  (<J)0/2)  be  much  larger 
than  unity. 

Comparing  Eq.  (C2)  with  Eq.  (Cl),  we  see  that  the  nonuniform  TE  current  predicted  by  the 
incremental  half-plane  strips  of  current  is  a  reasonable  approximation  to  the  actual  TE  current  when 
ld/ql  >  1/4;  that  is.  when  the  distance  from  the  leading  edge  along  the  diffracted  ray  is  less  than  about 
one-quarter  of  the  focal  length.  Thus,  secondary  diffraction  at  the  trailing  edge  of  the  nonuniform 
current  emanating  from  the  leading  edge  will  be  accurately  determined  by  the  half-plane  incremental 
current  strips  when  the  dislance  along  the  diffracted  ray  is  less  than  about  one-quarter  of  the  focal 
length.  Of  course,  for  a  plane  wave  incident  along  a  direction  far  from  grazing,  the  secondary 
diffraction  of  the  nonunlform  current  will  become  negligible  with  increasing  plate  size.  Moreover,  as 
mentioned  above,  near  grazing  Incidence  the  incremental  nonuniform  current  closely  approximates 
the  actual  current  across  the  plate,  and  secondary  diffraction  Is  accurately  predicted.  Interestingly, 
near  grazing  incidence  the  focal  length  of  the  diffracted  rays  approaches  infinity.  Hence,  although  the 
expressions  for  current  in  Eqs.  (Cl)  and  (C2)  are  not  valid  close  to  grazing,  they  correctly  predict 
equality  at  grazing  incidence  between  the  actual  and  incremental  half-plane  strip  nonuniform 
current. 
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